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Abstract 

Mathematical ecology and mathematical epidemiology are major fields 
in both biology and applied mathematics. In the present paper, a four- 
dimensional eco-epidemiological model with infection in both prey and preda- 
tor populations is studied. It consists of susceptible prey, infected prey, sus- 
ceptible predator, and infected predator. The functional response is assumed 
to be of Lotka—Volterra type. The behavior of the system such as the ex- 
istence, boundedness, and stability for solutions and equilibria are studied 
and also the basic reproduction number for the proposed model is computed. 
Moreover, a related control model and optimal treatment for the control 
model are presented. Finally, to verify the analytical discussion, a numerical 
simulation is carried out. 
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1 Introduction 


Health is the most important topic in the world for human, animal, plant, and 
any animate entity. Anybody gets involved with many various viruses such 
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as Black Death (1300), Spanish Flu (1918), HIV/AIDS (1920), Swine Flu 
(2009), Ebola virus (2013), Zika virus (2016), Corona viruses such as: SARS- 
Cov(2003), MERS-Cov(2012), COVID-19 (2019), and so on, which are the 
most deadly and disastrous viruses in 20th and 21st centuries. Some of them 
are not only epidemic but also are pandemic (that is a world epidemic). 
Black Death is just the most famous epidemic historically. The history of 
epidemic is an ever fascinating area. Hence nowadays, the study of diseases 
and their transmission is most studied subject. Therefore, physicians have 
to apply mathematician knowledge to protect. man against infected diseases 
transmitted by viruses. It is clear that ecology, biology, and epidemiology are 
major areas. Therefore, many researchers join together to work on common 
topics and interdisciplinary subjects. The above sentences motivate mathe- 
maticians, physicians, researchers, and scientists to join together in the said 
areas such as dynamics of infectious disease by the help of theory of modeling, 
dynamical systems, and so on. 

Infectious diseases may be a major factor in human population size. One 
of them is Black Death in the 14th century, which killed up to one-fourth of 
the people (which occurred in Eurasia); see [6]. Diseases of European people 
such as smallpox that were carried out by Cortez and coworkers to Mexico, 
decimated native population there in the 16th century. Rinderpest caused 
high mortality in wild animals in Africa at the end of the 19th century. 

It is obvious that describing any type of ecological phenomena such as 
predator-prey, coexistence, populations reactions, and the like entails having 
sufficient mathematical knowledge. Mathematical models may be applied 
as an effective tool to analyze such realistic phenomena. In the complex 
ecosystem some reactions such as coexistence, competition, food chain, food 
web, and predator-prey relationships may be important in regulating numbers 
of population of species. Whenever a bounty was placed on natural predators 
in the Kaibab Plateau in Arizona, the number of deer population went up for 
food supply, and then over half of the deer died in 1923-25; see [6]. Regarding 
to new researches, we may see an increase in a number of works that describe 
the relationships between demographic process among different populations 
and disease. Mathematical biologists indeed have been working on merging 
ecology and epidemiology (see [2, 7, 8, 5, 10, 12, 14, 15, 16]). 

This research work describes a predator and prey type by using biological 
and mathematical modeling. This model consists of two species: prey species 
and the predator species. Any species has two subclasses: susceptible and 
infected. 

The original parts of the paper are constructed as follows: Section 2 
consists the mathematical modeling; indeed a model is formulated. Section 
3 involves stability analysis consisting boundedness, equilibria, stability, and 
the basic reproduction number Ro for model (2). In section 4, we analyze 
a control model having infection. Also an optimal treatment approach is 
studied. Indeed, this section consists making control model for system (2), 
applying the control functions for system (3) and its analysis, and establishing 
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an optimal control to minimize the number of infected species. In section 5, to 
verify the analytical discussion, numerical simulations are computed. Finally, 
in the last section conclusion is presented. 


2 Mathematical modeling 


In this section, a predator-prey system having infection is described and then 
its mathematical model is presented; see [1]. First assume that S(T), I(T), 
X(T), and Y(T) describe the densities of susceptible prey, infected prey, 
susceptible predator, and infected predator at time 7’, respectively. In the 
second step, some basic hypotheses which have been made in the formulating 
of model, are presented as follows: 

H,) The prey species grow logistically having intrinsic growth rate y and 
carrying capacity K (which is positive) in the case when there is no diseases 
and predator species. 

Hy) Just susceptible prey may reproduce. 

H3) Prey species can be infected never recovers, which will either die or be 
removed in predation interacting. The infected prey species have a disease 
induced death. 

H4) Infected predators species can catch a healthy predator. Meanwhile, 
infected prey may be available for predation by predators. 

H;) By considering a direct contact with infected predator, disease may be 
spread. 

Hg) Finally, consider predator species die naturally, but infected species have 
a disease induce excess death rate also. 


These mentioned considerations motivate one to the following system: 


dS S+I 


= 1 T—b,SX 

dT 75( K )= a8. 1SX, 

os Sf —dI —filxX —mlY. 

dT ay 1 1 MyLr, 

dX 

ar — cy SX } glx eyXY 6, X, 

dY 

ar =e XY — (d1 +ayz)Y + n11Y, (1) 


with the following initial values S(0) = SoI(0) = Ip, X(0) = Xo, Y(0) = Yo. 
We moreover make an obvious assumption that all the parameters are 
positive. Indeed, the parameters are described as follows: 


7) a, and 6, are infection rates for prey population and predation rate of 
susceptible prey by healthy predators. 


it) d,is an infected prey species disease induced death. 
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iti) f; and m, describe predation rate of infected prey by susceptible and 
infected predators, respectively. 

iv) c, and g; describe conversion rates for healthy and infected preys to 
healthy predator, respectively. 

v) 01 describes the natural death rate for the predator population having a 
disease induced death rate a}. 

vt) e; and n; denote infection rates for predator population and conversion 
rate for infected prey to infected predator. 

By using the following scaling 

8 3 yt 4a x iY x, and t = yT, after some simplification, system 
(1) reads as follows: 


d 

G7 8k - 8-8) — asi — bse, 

di 

a = asi — di— fxi— miy, 

dx _ 4 gi 5 

qe 7 Coe t git — exy — 5a, 

d 

7 =exy —ayt niy. (2) 


_ ak _ b1K _ ak _d _ ek — fik 
arr ear CS a d=, ee f= y? 
K 9 
g= gi , m= malt n= mk a= Ste. i= cm 
a x 9 Y Y 


We are now going to study the solutions boundedness for system (2). For 
system (2), any solution that starts in R4, remains in this area. 

As regarding any parameters nonnegative, one is able to see that the 
right-side of model (2) can be a smooth function at (s,i,2,y) in the first 
octant, 

NX = {(s,i,2,y)|s > 0,i > 0,2 > 0,y > O}. 


Also Q is an invariant set. 

Since system (2) is homogeneous, we have s = 0,7 = 0,2 = 0, and y = Ois 
the trivial solution. The uniqueness and existence theorem ensures that any 
trajectory starting from the first quadrant remains in it, that is, no trajectory 
will cross the coordinate planes. 


3 Stability analysis 


In the present section, the solution behavior such as boundedness and sta- 
bility for equilibria, is studied. Then, the basic reproduction number for the 
said model is computed. 
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3.1 Equilibria Stability 


Theorem 1. [1] Any solution for system (2) is uniformly bounded. 

In the present section, first the equilibrium points and also the basic 
reproduction number Rp are computed for system (2). Then the stability 
analysis for these equilibria will be presented. Now, we are going to find out 
equilibria for system (2) as follows: 

a) Eo(0,0,0,0), which is known as the trivial equilibrium. 

b) £,(1,0,0,0), which is known as the axial equilibrium. 

c) E2(s2,0,x%2,0), which is known as the disease-free equilibria. In the 


truth of the following conditions, this equilibrium exists: c > 6, sz = 7 and 
m= ae Moreover, this condition c > 6 indicates that c,; > 7 that is, 


a disease-free equilibrium exists provided that predator is a high capacity 
consumer. 

d) £3(s3,%3,0,0), which is known as the predator-free equilibria; 
In the ee of the following conditions, this equilibrium exists a > d, s3 = d 


, 13 = a Furthermore, this condition a > d implies a, K > dj, that is, 
it indicates that the predator-free equilibrium exists provided that a disease 


induced death rate is low enough. 


e) E4(84, 14, £4, 0), which is known as the infected-predator-free equilibria. 
This exists when s4,%4,2%4 are given by 


g(bd + f) — iti ) 
g(ab+ f)—e(f +af)’ 
_ d(ab+ f) — c(bd + f) 
glab+ f)—e(f +af)’ 


_ ag(bd + f) — gd(ab+ f) + (f +af)(cd — ad) 
gf(ab+ f) —cf(f +af) 
f) E5(s5,0,25, ys), which is known as the infected-prey-free equilibria. 


In the case of e > ba and ce > bca + ed, this equilibrium exists, where 
$5,25,Y5 are given as 


e — ba 
685 = ) 
€ 
a 
w=; 
e 
ce — bca — ed 
rr 
e2 


g) E*(s*,i*,x2*,y*), which is known as the interior equilibrium point; it 
exists provided that 
(i) |Po| > max{|(1 + a)P, + bPg|, tee oFr thers |) 
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(ii) P,, P2, and P3 are of the same sign, 


where P,, P2, and P3 are given by 


P, = e(aba + fa+em + de) — (bema + med + ae’), 


Py = e(abn + fn + (1+ a)em) — (bemn + a(1 + a)e” + meg), 
P3 


(aen + mné + (1+ a)cema) — cmn + den + (1+ a)aea + mga). 


If the above conditions are true, then the values of s*,7*, x*, y* may be found 
by 


 — faa +a)P — bPs 


P, 
a 1 
P,’ 
Pees 
= 5 
» _ (g—c—ac)P, + (c— 6)P2 — beP3 


eP» 


Now we are able to analyze the stability of the above equilibrium points 
for system (2). Indeed, some conditions are presented which show the points 
Ey, Fo, E3, E4, Es and E* are locally asymptotically stable. 


Theorem 2. [1] The following statements are true for system (2). 
(a) The trivial equilibrium Ep is unstable. 


(b) If d > a and 6 > cc, then the equilibrium Fj is locally asymptotically 
stable. 


(c) If e6+bca > ec, A, > 0,A3 > 0, and A; Ag > As, then the equilibrium 
Ey is locally asymptotically stable. 


(d) If n(a — d) < aa(1 + a)andg(a — d) + (cd — ad)(1 + a) < 0, then the 
equilibrium EF3 is locally asymptotically stable. 


(e) Ifc < 6, B, > 0, Bs > 0, and B, By — Bs > 0, then the equilibrium EF, 
is locally asymptotically stable. 


(f) If Cy > 0,C3 > 0, C4 > 0 and CyC2C3 > ce + CPC, then the 
equilibrium Es is locally asymptotically stable. 


(g) If D; > 0,D3 > 0,D4 > 0 and D,D2D3 > De as D?Da,, then the 
equilibrium E”* is locally asymptotically stable. 
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3.2 The basic reproduction number 
The next generation matrix method [3] is used to calculate the basic re- 


production number Ro; see [9]. Clearly, i and y are the relevant classes of 
infection. The classes i(¢), y(t) from our model are 


di 
im di — fai — miy, 
dy =exy—ayt+niy 
dt ae 


To define the basic reproduction number, denote the vectors y and v and 
the inflow and outflow from the disease compartment 7 and y as follows: 


_f{gi\_ ast 
 \ pe) ~ Vey +niy )? 
ec PE Ve di+ fai+myt 
le J ay : 


We compute the 2 x 2 Jacobian matrices, evaluated at the disease- free 
equilibrium point £2(s2,0, x2, 0), 


F = ( $24(s2,0,2r2,0) ) 


and 


2 


V = (Bet (s2,0,2r2,0) ) 


Therefore, two matrices F’ and V corresponding to the gain and loss 
components of system (2) are defined as 


as 0 asz O 
F= ; = 
ny ex + ne ” 0 ex, 
2 


ve Gare 
|B. 


and 


0 a 0 a 


Now, the next generation matrix is defined as 


as2 0 
— ji Lo 
G=FV'= ( bar 2) ; 


(exgd + ex? f + asga) €Xx2089 
a(d + fx) 


2 
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The basic reproduction number indeed may be dominant eigenvalue for the 
next generation matrix. Thus, 


R= 1 (exod + exo? f + asea) _ | (ease)? + (ex2)*(d + fae)? 4 ©%2082 
aa a(d + fx) | a?(d+ fae)? a(d + fxg)”’ 
where 
; ) : c—6 
2 c’ 2 > be 


If Ro > 1, then the disease is endemic. If Rp < 1, then the disease-free 
equilibrium is locally asymptotically stable. 


4 Control model and optimal treatment 


In this section, first we make a control model for system (2) and then analyze 
the application of control for system (3). After it we describe an optimal 
control approach. 


By this mean, we add the control functions uw; and uz as follows: 
(i) uz is denoted as the control factor for the first and second species (prey 
and infected prey species). 
(ii) uz is denoted as the control factor for the third and fourth species (preda- 
tor and infected predator species). 
The role of control functions is to treat the infected prey and predator species. 
Hence, the control system can be written as follows: 


ds 

dt = s(1—s—1) —a(1—u))si — bsz, 

di 

a = all u1)si — di — fai — miy, 

dx 

qe OSE + giz — e(1 — ug)ry — da, 

d 

7 = e(1 — ug)xy — ay + niy, (3) 


where 0 < u; < 1,i = 1,2. When u; = 0, no treatment occurs and when- 
ever u; = 1, the model shows the full treatment. To investigate the model, 
among all of equilibrium points, we analyze the interior equilibrium, that is, 
E*(s*,i*,x*,y*) for system (3). It is easy to show that E* exists, if 


(i) (Py) > max{|(1 + o(1 — a) P; + BP,|, ee ere teers 


(ii) Pi, P2, and P3 are of the same sign, where P;, P:, and P3 are given 
by 
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P, =e(1 — u2)(a(1 — u1)ba + fa + com + de(1 — ug)) 
— (bema + me(1 — ug)d + a(1 — uy)e?(1 — u2)?), 

Py =e(1 — ug) (a(1 — u1)bn + fr + (1 4+ a(1 — ui )em) 
— (bemn + a(1 — u)(1 + a(1 — uz))e?(1 — ug)? + me(1 — u2)g), 

P3 =(a(1 — ui )e(1 — ug)n + mnéd + (1 + a(1 — u1))ema) — emn + de(1 — u2) 
+ (14 a(1 — u1))a(1 — ur )e(1 — u2)a + mga). 


Ps) 


When the above equalities are satisfied, the values of s*,7*,x*,y* are as 
follows: 


s* Pp — (1+ a(1—u1))P, — P3b 


P» : 
settee 
ae. 
18 
v= Py’ 
«  (g-—c-—a(l—wui)c)P; + (c — 6)P2 — bcP3 
= e(1 — ug) P» ; 


Theorem 3. The nontrivial equilibrium point E* is locally asymptotically 
stable for model (3) provided Di > 0,D5 > 0,D, > 0 and D{D,D5 > 
Di DI Di. 


Proof. At the interior equilibrium E*, the Jacobian matrix J(£*) can be 
obtained as follows: 


—s* —(1+a(1—1u))s*  —bs* 0 
x, | a(1—w)i* 0 —fa* —mi* 
JB") = ca* gu* 0 —e(1 — u2)a* 
0 ny* e(1 — u2)y* 0 


The corresponding characteristic equation is given by 


do DLA + DEA? + DIA+ Di =0, 


Di =e?(1 — up)?x*y* + fgx* + a(1 —u)[a(1 — uw) + 1]s*i* + mni*y* 
+ bes*x*, 

Di, =e?(1 — ug)2s*a*y* + fgs*x*? + mns*i*y* — e(1 — u2)nfa*?y* 
+me(1 — u2)gi*x*y* + a(1 — u)bgs*i*a* — cf(a(1 — uy) + 1)s*2"*, 


and 
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Di, =—e(1 — ug)nfs*x**y* + [e(1 — u2)mg 
+ a(1 — uz) (a(1 — uz) + 1)e(1 — ue) — a(1 — ui) be(1 — u2)n 


OK OK ok 


+ bemn — cme(1 — u2)[a(1 — ui) + 1)]s*i*a*y*. 


Regarding the Routh—Hurwitz criterion, it is easy to see that if 


Di>0, i=1,3,4, 
Di (D),D} = D,D) —D? > 0, 


then all of the eigenvalues for J(E*) have negative real parts. 
Therefore, if D, > 0, D4, > 0, Di, > 0 and D/,D4,D!, > D4? + Di? D4, then 
E* is locally asymptotically stable. 


Now, we describe an optimal control approach. As it is well known, the 
optimal control problem consists of forcing the solution for a system out of 
undesirable set at a given time. We are going to control the population of 
infected prey and predators, in an optimal method. And so, we establish an 
optimal criterion that consists of the minimizing of total number for infected 
species along with the treatment costs. We need to find the control functions 
U1, U2, as defined before, such that minimize the objective functional 


T 
min f [5 (Wien? + Woug”) + i(t) + y(t)]dt. (4) 


The existence of an optimal control for the state system is analyzed by using 
the theory developed by Fleming et al. [?]. The boundedness of the solu- 
tions was discussed in section 3, which is needed to obtain the existence of 
an optimal control. Here, we may state the existence theorem as follows: 


Theorem 4. There exist optimal controls u; and ug that minimize the 
objective functional (4), if the following conditions are met: 


1. The class of all initial conditions with controls u and v such that u 
and vu are Lebesgue integrable functions on [0,7] with values in the 
admissible control set along with each state equation being satisfied is 
not empty. 


2. The admissible control set is closed and convex. 


3. The right-hand side of the state system is continuous, bounded above 
by a sum of the bounded control and the state, and it can be written 
as a linear function of u, and ug with coefficients depending on time 
and the state variables. 


4. The integrand of the functional is convex on the admissible control set 
and is bounded above. 
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The classical approach is described in [11], so we can derive the associated 
Hamiltonian as: 


H = 5 (Wu? + Woug”) + i(t) + y(t) + ps(t)(s(1 — s — i) — a(1 — uw) si — bsz) 


+ p;(t)(a(1 — ui) si — di — fai — miy 
+ px (t)(csa + gia — e(1 — ug)xy — 62) 


( 
+ py(t)(e(1 — ua)ay — ay + niy). 


where the functions p(t), p(t), px (), py(t) are called as the co-state variables. 


These variables should satisfy in the following equations: 
dp,(t) OH  pj(t) OH dp,(t) OH dp,(t) OH 


dt Os’ dt Or”? dt Ox’ dt Oy’ 
or explicitly 


dp.(t) OH 
dt Os 
= —(p,(t) — 2sp,(t) — ip,(t) — aisp,(t) + aisuip.(t) — bxp,(t)), 
ait) - see = —(1— sp,(t) — asp,(t) + asurps(t) + asp;(t) 


— asuyp;(t) — dp;(t) — fapi(t) — mip;(t) + gxp.(t) + nypy(t)), 


Walt) _ _ OM _ _(_bsp,(t) — fipslt) + espe(t) + gipe(t) ~ evpelt) 
+ eyurpe(t) — dpe(t) + eypy(t) — uzeypy(t)), 
tpult) = a = —(1 — mip;(t) — exp,(t) + exuepz(t) + expy(t) 
— exUuzpy(t) — apy(t) + nipy(t)). (5) 


For the above system of co-stated variables (5), we consider the following 
boundary conditions: 


s(0)=5., i(0)=%, 2«#(0)=2., y(0) = Yy%, 


and 


ps(T) = pi(T) = p2(T) = py(T) = 0. 


Differentiating H with respect to wu, and ug yields: 


= = Wiu + asi(ps(t) — pi(t)) = Wius + (asi)ps.(t) — (asi)p,(t)) = 0, 
Fa = Wou2 + exi(pz(t) — py(t)) = Waue2 + (ext)p,(t) — (ext)py(t)) = 0. 


Thus, the control factors u; and wz can be obtained as follows: 
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u(t) = F(-palt) + vilt) 
ualt) = (Pelt) + py(t)) 


As regards 0 < uj, ug < 1, we have the optimal control laws as 


u,*(t) = min{1, max{0, ui}} = min{1, max{0, Wt) — p,(t))}}, 
ug*(t) = min{1, max{0, u2}} = min{1, max{0, Full — px(t))}}. 


5 Discussion and numerical simulation 


It is obviously that the numerical simulation is important beside the analyt- 
ical findings to verify them. By using MATLAB software, we now present 
computer simulation for different solutions of systems (2) and (3). We first 
present a computer simulation for equilibrium of system (2). The behavior 
of these equilibria are shown in Figures 1-6. By setting the parameters and 
initial conditions for system (2), we can construct Tables 1 and 2, which are 
referred in Figures 1-6. These values are calculated in MATLAB software. 

Now, we are going to illustrate Figures 1-6. First, we take the param- 
eters and initial condition for system (2) as Row No. i of Tables 1 and 2. 
Then we sce that the conditions of Theorem 2 are satisfied and consequently 
£,(1, 0,0, 0) is locally asymptotic stable. This illustration is shown in Figure 
la. After it, we consider parameters for system (2) as Row No. ii of Table 1 
and Row No. 1 of Table 2, then the conditions of Theorem 2 are satisfied and 
consequently F'2(0.8,0,0.2,0) is locally asymptotic stable. The behavior of 
5,1, x,y with t is depicted in Figure 1b. Then, we take parameters and initial 
condition for system (2) as Row No. ii) in Table 1 and Row No. i Table 
2, then the conditions of theorem 2 are held, and so E3(0.25, 0.5357, 0,0) is 
locally asymptotic stable. The behavior of s,i,x,y with t is depicted. This 
description is plotted in Figure 1c. In the next step, we take parameters and 
initial condition for system (2) as Row No. iv in Table 1 and Row No. ii of 
Table 2. 

Therefore, the conditions of Theorem 2 are satisfied, and so 
4 (0.4884, 0.2519, 0.1589,0) is locally asymptotic stable. This locally be- 
havior is shown in Figure 1d. Then, If we take the parameters and initial 
condition for system (2) as Row No. v in Table 1 and Row No. itz of Table 
2, then the conditions of Theorem 2 are satisfied and therefore the equilib- 
rium point £5(0.75,0,0.25, 0.4583) is locally asymptotic stable. This kind 
of stability is plotted in Figure le. Finally, we consider parameters as Row 
No. vi in Table 1 and Row No. iv of Table 2. Then the equilibrium point 
E5(0.75, 0, 0.25, 0.4583) is locally asymptotically stable and s,i, x, y approach 
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to (s*,i*,x2*,y*) = (0.2568, 0.1921, 0.1550, 0.1713), which are shown in Figure 
iy. 

Now, we can solve our optimal control problem by the help of an iter- 
ative method numerically. We here obtain the optimality system from the 
state and adjoin equations. The optimal control problem strategy is obtained 
by solving the optimal system, which consists of eight ordinary differential 
equations and boundary conditions. By the Runge-Kutta fourth scheme, 
we can solve the optimality system; see [4]. Starting with an initial guess 
for the adjoin variables, the state equations are solved by a forward Runge— 
Kutta fourth scheme in time. Thus, by a backward Runge-Kutta fourth 
scheme, those state values are used to solve the adjoin equations because 
of the transversally conditions; see [13]. We present the results using the 
"bvp4c” and "bvp5c” subroutines in MATLAB software. To control in all in- 
fected prey and predator individuals, we use two control factors as treatment. 
We consider the treatment for 30 days, because a long treatment in the form 
of medication has a potentially harmful side effect and the best time of vacci- 
nation is the possible early stage of diseases. We here use a set of parameters 
value a,b,c,d,e, f,g,m,n,a,6 to determine the numerical simulation of the 
optimality system with a sufficiently small time step size. In Figures 2-7, we 
plot all population sizes in two systems having control and without control. 
The solid line denotes the population of individuals in the system without 
control while the dash-dotted line denotes the individuals population in the 
system with control. The changes of population is clearly visible for equi- 
libriums F3, E,, and E* in Figures 4, 5, and 7. The population of infected 
individuals, i(t) and y(t) is reducing during the time period, in these cases. 
Moreover, the population of the other susceptible individuals, s(t) and x(t) 
is increasing. 

The equilibria EF, and £2 are not infected, then the control factors have 
no important effect on results as we see in Figures 2 and 3. In the other 
hand, after applying the control function F3, s3(t) increases significantly 
and 2(t) decreases. Figure 3 demonstrates the increase of prey and decrease 
of infection in prey population. Also, it reveals less decrease in predator 
population. After applying controls, Figure 7 at the E* equilibrium shows 
the decrease at infected prey and predator, while increasing the predator 
populations. 
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(€) Row No. v of Table 1 and Row No. iii of Table 2 (£) Row No. vi of Table 1 and Row No. iv of Table 2 


Figure 1: Stability of system (2) for different parameters and initial conditions. 
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Table 1: Parameter values for system (2) 


Row 
No.Parameter a b c d € f g min a 6 
i 0.4} 1 | 0.1) 0.5 | 0.2 | 0.2 | 0.15 | 1.5) 1.4 | 0.5 | 0.2 
ii 04/} 1 | 0.5 | 0.5 | 0.2 | 0.2 | 0.15 | 1.5 | 1.4 | 0.5 | 0.4 
iii 0.4} 1 | O.1 | OL | 0.2 | 0.2) 0.15 | 1.5 | 1.4 | 0.8 | 0.4 
iv 04/ 1 | 01} 01) 02)06)] 06 }15),14) 1 | 0.2 
v 0.4} 1 1 | 01 )1.2/06/] 04 | 1.5] 1.4) 0.3 | 0.2 
vi 1.9 | 1.2) 0.8 | 0.2 | 0.2 | 0.2 | 0.15 | 1.5 | 1.4 | 0.3 | 0.2 


Table 2: Initial conditions for system (2) 


Row 
No, Initial value s(0) i(0) x(0) y(0) 
i 0.5 | 0.5 | 0.5 0.5 
ii 0.5 | 0.5 | 0.5 0.2 
ill 0.5 | 0.2 | 0.5 0.2 
iv 0.5 | 0.2 | 0.5 0.5 


Populations 


Figure 2: This plot shows the population of s(t), i(t), x(t), y(¢) for equilibrium EF) with 
control (solid line) and without control (dash-dotted line). 
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Populations 


Figure 3: This plot shows the population of s(t), i(t), x(t), y(t) for equilibrium 22 with 
control (solid line) and without control (dash-dotted line). 


Populations 


Figure 4: This plot shows the population of s(t), i(t), x(t), y(t) for equilibrium £3 with 
control (solid line) and without control (dash-dotted line). 
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Populations 


Figure 5: This plot shows the population of s(t), i(t), x(t), y(t) for equilibrium 24 with 
control (solid line) and without control (dash-dotted line). 


Populations 


Figure 6: This plot shows the population of s(t), i(t), x(t), y(t) for equilibrium 25 with 
control (solid line) and without control (dash-dotted line). 
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Figure 7: This plot shows the population of s(t), 7(t), x(t), y(t) for equilibrium E* with 
control (solid line) and without control (dash-dotted line). 


6 Conclusion 


Last century was known for mathematical physics era, but the present century 
belongs to biological mathematics, which has achieved considerably advanced 
developments in all branches of bioscience, medical, ecology, and so on. In the 
present research work, a four-dimensional Eco-epidemiological model having 
infection in the present species is studied. In other words, we formulated a 
prey-predator model having a disease in both of the populations. For almost 
all the models with diseased prey, we considered the predators live only on 
the infected preys (as they are weak and more vulnerable), therefore, the 
susceptible preys are completely out of danger is an oversimplification. In 
our model, we made more realistic assumptions which are as follows: 


(i) The susceptible predators are capable of catching both the susceptible 
and infected preys, and 


(ii) infected predators (being weak with disturbed internal mechanism) 
can manage only the infected preys (due to the same reasons). 


The model may be used in many ways. Precisely, the behavior of solutions 
and equilibria such as existence, boundedness, and stability are investigated. 
Meanwhile the basic reproduction number for the proposed model is com- 
puted, the related control model and an optimal treatment are carried out. 
New controlled model are developed from the numerical simulation of the 
optimal system, which represents dynamics in each individual of the commu- 
nity. The controls developed in this paper support that the number of the 
susceptible and infected individuals decreases and the number of the recov- 
ered individuals increases in the optimal system. We also pointed out that 
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for certain values of the control rate, there exists its corresponding optimal 
solution. The above discussion is verified by the numerical simulation. 


Acknowledgment 


Authors are grateful to there anonymous referees and editor for their con- 
structive comments. 


References 


10. 


. Bera, S.P., Maiti, A., and Samanta, G.P. A prey-predator model with 


infection in both prey and predator, Filomat 29(8) (2015), 1753-1767. 


. Cojocaru, M.G., Migot, T., and Jaber, A. Controlling infection in 


predator-prey systems with transmission dynamics, Infect. Dis. Model. 
5 (2020), 1-11. 


. Diekmann, O., Heesterbeek, J. A. P., and Roberts, M. G. The construc- 


tion of next-generation matrices for compartmental epidemic models, J. 
R. Soc. Interface, 7 (2010), 873-885. 


. Heath, M.T. Scientific computing: An introductory survey, The McGraw- 


Hill, New York, 2002. 


. Hethcote, H.W. The mathematics of infectious diseases. SIAM Review, 


42 (2000), 599-653. 


. Hethcote, H.W., Wang, W., Han, L., and Ma, Z. A predator-prey model 


with infected prey, Theor. Popul. Biol. 66 (2004), 259-268. 


. Jana, S. and Kar, T.K. Modeling and analysis of a prey-predator system 


with disease in the prey, Chaos. Soliton. Fract. 47 (2013), 42-53. 


. Jiao, J.J., Chen, L.S., Nieto, J.J., and Angela, T. Permanence and 


global attractively of the stage-structured predator-prey model with contin- 
uous harvesting on predator and impulsive stocking on prey, Appl. Math. 
Mech-Engl. 4(3) (2008), 181-191. 


. Jones, J.H. Notes on Ro, Standford University, 2007. 


Hugo, A. and Simanjilo, E. Analysis of an eco-epidemiological model un- 
der optimal control measures for infected prey, Appl. Appl. Math. 14(1) 
(2019) ,117-138. 


176 


11. 


12. 


13. 


14. 


15. 


16. 


M.H. Rahmani Doust, M. Shirazian and M. Shamsabadi 


Kamien, M.I. and Schwartz. N.L. Dynamic optimization: The calculus of 
variations and optimal control in economics and management. Courier 
Corporation, 2012. 


Kumar, S. and Kharbabda, H. Stability analysis of prey-predator model 
with infection, migration and vaccination in prey, arXiv:1709.10319 
(2017). 


Lenhart, S. and Workman, J.T., Optimal control applied to biological 
models, Chapman and Hall/CRC, Mathematical and Computational Bi- 
ology Series, London UK, 2007. 


Pan, S. Asymptotic spreading in a lotka-Volterra predator-prey system, 
J. Math. Anal. Appl. 407(2) (2013), 230-236. 


Simon, J.S.H. and Rabago, J.F.T., optimal Control for a predator Prey 
Model with disease in The Prey Population, Malaysian J. Math. Sci. 12(2) 
(2018), 269-285. 


Zhou, J. and Shi, J. The existence, bifurcation and stability of positive 
solutions of a diffusive Leslie-Gower predator-prey model with Holling- 
type II functional responses, J. Math. Anal. Appl. 405(2), (2013), 618- 
630. 


